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<— ! Abstract 

We consider a semiconductor in a non-equilibrium steady state, 

with a dc current. On top of the stationary carrier motion there are 

fluctuations. It is shown that the stationary motion of the carriers 

(i.e., their drift) can have a profound effect on the electromagnetic field 

fluctuations in the bulk of the sample as well as outside it, close to the 

Q surface (evanescent waves in the near field). The effect is particularly 

pronounced near the plasma frequency. This is because drift leads to 

a significant modification of the dispersion relation for the bulk and 

r"^ surface plasmons. 

> 

^ 1 Introduction 

CN 

Random thermal motion of charge carriers in a body produces a fluctuating 
O electromagnetic field. Properties of such a fluctuating field have been stud- 

ied for a very long time and are discussed in a number of textbooks [1-4]. 
Outside the body one should distinguish between the near and far field do- 
^ main. In far field, i.e., when the distance £ from the surface is much larger 

than the wavelength A of the corresponding Fourier component of the field, 
one observes the well known phenomenon of thermal radiation. In the oppo- 
site regime, I << A, there exists (in addition to radiation) a non-radiative 
electromagnetic field that is due to evanescent waves excited by the jiggling 
carriers. This random evanescent field, close to the sample surface, is ro- 
tationless and can be described by a scalar potential. Thermal fluctuating 
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fields manifest themselves in a variety of experimentally observable physical 
phenomena (the Casimir-Lifshitz forces, near-field heat transfer, noncontact 
friction) which motivates the ongoing interest in the subject (see [5-8] for 
recent reviews). 

In the present paper we consider an out of equilibrium, steady state situ- 
ation when a dc current is established in a conducting medium. Drift of the 
carriers can have a profound effect on the fluctuating electromagnetic field 
inside, as well as outside, the medium. There exists a large body of work on 
current fluctuations in the presence of drift (for some early references see [9- 
13]). There is clearly a connection between that old work and the subject of 
this paper, in which we emphasize some conceptually important points and, 
in particular, consider the effect of drift on the fluctuating evanescent field 
outside the medium, close to its surface. It is known that drift of the carriers 
affects the field outside the medium and can even lead to the existence of a 
new type of weakly decaying surface wave [14,15]. The possible effect of such 
waves on the electromagnetic field fluctuations close to the sample surface 
was pointed out long ago [16]. Here we take a broader view of the problem 
and consider a more general model of a conducting medium. This enables us 
to discuss the fluctuational field at frequencies higher than the inverse scat- 
tering time of the carriers, when plasmonic excitations come into play. Drift 
is expected to have a particularly strong effect near the plasmon frequency. 

The organization of the paper is as follows: In Sec. 2 we formulate the 
model and derive the corresponding permittivity tensor, with respect to the 
steady state. This permittivity tensor relates quantities which fluctuate on 
the background of the stationary motion of drifting carriers. In Sec. 3 we 
summarize Rytov's method for treating the electromagnetic field fluctuations 
and introduce modifications needed for application of the method to our prob- 
lem. In Sec. 4 we study the effect of carrier drift on the field fluctuations well 
inside the sample (the limit of an infinite medium). The experimentally more 
relevant case of a field outside the sample, close to its surface, is considered 
in Sec. 5. It is shown there that drift of the carriers leads to a pronounced 
dip in the field power spectrum in the vicinity of the bulk plasma frequency. 
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2 Permittivity tensor and plasma waves in the 
presence of drift 



We consider a conducting medium, e.g., a semiconductor, subject to a con- 
stant electric field E . This field causes drift of the carriers, with the charge 
e, so that there is a steady state current density j = enov , where n is the 
equilibrium density of carriers and v is their drift velocity. On top of the 
stationary motion there are fluctuations. All fluctuating quantities will be 
denoted by the corresponding letters without any subscript or superscript. 
For instance, E(r, t) and j(r, t) represent fluctuating parts of the electric field 
and current density at point r and time t. Relations between the fluctuating 
parts of various quantities are obtained by linearization near the steady state. 
Since we shall be dealing with temporary and spatial dispersion, algebraic 
relations (rather than integral ones) exist only for the Fourier components. 
In particular, 



where tilde indicates the Fourier transformed quantities and cr a p(oj, k) is the 
conductivity tensor with respect to the steady state (summation over j3 is 
implied). The specific form of this tensor depends, of course, on the model 
used for the carrier motion (see below). 

Let us emphasize that current density j accounts only for the motion of 
the mobile carriers but not for the polarization current of the lattice. The 
latter is incorporated into the dielectric constant, e^, of the lattice, so that 
the electric displacement (its fluctuating part) is 



where e a p(u), k) defines the steady state permittivity tensor. We assume that 
the frequency oj is far from any resonant frequencies of the lattice and set 
€l = const, thus neglecting any possible dispersion effects in the lattice. 

We are interested in the effect of the longitudinal plasma waves on fluctua- 
tions. Therefore we consider a rotationless electric field, E(r, t) = — V$(r, t), 
which in the absence of sources satisfies 



(1) 



D a (uj, k) = e L E a (u, k) + i — a a/3 (uj, k)Ep(u, k) 



UJ 



e al3 (uj,k)Ep(uj,k), (2) 




(3) 
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where the caret emphasizes that, in the presence of dispersion, e a p is an 
integral operator relating the electric displacement at point r and time t to 
the electric field at earlier times in some vicinity of r. Fourier transforming 
(3) yields 

k 2 e(ou,k)$(uj,k) = 0, (4) 

where a scalar quantity 

e(w,k) = ^^e a(3 (u,k) (5) 
has been introduced. The equation 



e(w,k) = (6) 

defines the dispersion relation for longitudinal waves in the medium [17]. 

In order to obtain an explicit expression for e Q/ g(u;,k) we will use a hy- 
drodynamic equation for the carrier flow V: 

+ ( v • V)V = -E - vV - -\v P , (7) 
at m mN 

where m is the effective mass of a carrier, z/V describes relaxation of velocity 
due to collisions, with frequency u, and the last term accounts for thermal 
pressure of carriers. Note that V,E and N refer to the total velocity, field 
and carrier concentration, i.e., V = v + v, E = E + E and N = no + n. 
The pressure is related to the concentration as p = NT, where T is the 
temperature in units of the Boltzmann constant ks- Thermal pressure does 
not play a major role in our considerations, which are focused on the effect 
of drift. It will be needed in the treatment of the infinite medium (Sec. 
5) because without thermal pressure (or some other mechanism of spatial 
dispersion, e. g. diffusion) one would encounter diverging integrals. Various 
versions of Eq. (7) are often used in semiconductor, as well as plasma, physics 
- see, e.g., Ref. 15 where the magnetic field effects are also included. Eq. (7) 
should be supplemented by the continuity equation 

dN 

— + div(iVV) = 0. (8) 
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Linearizing (7), (8), as well as the total current density eiVV, with respect 
to the fluctuating quantities n, v, E, and Fourier transforming to u>, k, one 
obtains 



j(w,k) = i 



.e 2 n 



m (3 + 



x 



X — {e(w, k) + (k • E( w , k)) 
T/m 



where 



/3 V ° + + -Tk 2 /m 



u — k • v . 



(9) 



(10) 



The conductivity tensor a Q( g(a;,k) is readily read off from (10). We write 
directly the permittivity tensor, as defined in (2): 



e aj g(w,k) = (e' L + ie'[)S a/: 



((3 + iv)u 



' 1 , Tk a k p /m 

<W + ^ + W 



IV) 



Tk 2 /m 



,(H) 



where w 2 = 47re 2 n /m and the lattice dielectric constant has been sep- 
arated into the real (e' L ) and imaginary (e" L ) parts. It follows from (5) and 
(11) that 



e(u;,k)=e'Jl 



w 2 



((3 + iv)u 



+ 



(3 13(13 + w) - k 2 R 2 D u 2 



(12) 



where co p = oo p / y e' L is the plasma frequency renormalized by the dielectric 
constant of the lattice, and R D = ^Jt/itlu 2 is the Debye screening radius. 
Neglecting the thermal pressure term one obtains 



e(u, k) = e' L 



1 



out. 



k • v + iv)(uu — k • v ^ 



(13) 



In equilibrium, i.e, for v = one recovers the standard Drude model (in the 
presence of the lattice): 



out 



(u + %v)uu e L 



(14) 



Since plasma waves will play an important role in our treatment of fluctu- 
ations, we pause to discuss briefly propagation of these waves in the presence 
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of drift. To see the effect of drift most clearly, let us neglect all dissipative 
terms [y — > 0, e" L — > 0) and compare the dispersion relation in equilibrium 
with that for vo 7^ 0. Using in (6) the expression (12), with v = e" L = v q = 
(i.e., (3 = u), we obtain the well known dispersion relation for the equilibrium 
plasma excitations: 

k 2 Rl = - 2 -l, (15) 

id 2 

which tells us that no wave can propagate for u < td p . For u > Co p propaga- 
tion is possible, provided that kRp is a small number - otherwise the phase 
velocity of the wave becomes close to the thermal velocity of the carriers and 
a strong collisionless (Landau) damping sets in [18]. Thus, the frequency of 
a propagating wave should be somewhat larger than the plasma frequency, 
id ~ Cu p (l + \k 2 R 2 D ). For the out of equilibrium situation, Vo 7^ 0, a similar 
calculation, under the conditions kRo < 1, k • v < u yields the dispersion 
relation 

id = u p + k • v + ^u p k 2 R 2 D , (16) 

so that, in contrast with the equilibrium case, propagation with frequen- 
cies below ujp becomes possible, if k • v is negative and sufficiently large 
in magnitude for (16) to be satisfied. The necessary condition for this is 
v > Rd^2uj p {£o p — uj). Although the effect of drift on plasma waves is 
an interesting topic in itself, we shall not pursue it any further but rather 
concentrate on the effect of drift on fluctuations. 



3 BASIC EQUATIONS OF THE THEORY 

In our treatment of fluctuations we follow Rytov's method [2-4], in which 
random Langevin sources are introduced into the Maxwell equations. These 
sources describe the spontaneous (thermal and quantum) fluctuations of po- 
larization and current density, j^(r,t). In equilibrium the correlation func- 
tion of these sources is determined by the fluctuation-dissipation theorem and 
it is given in terms of the imaginary part of the dielectric constant. Since 
we are dealing with a nonequilibrium situation, there is no general relation 
between the correlation function of the sources and either e a p((d,k) or its 
equilibrium counterpart. However, it can happen that, in spite of the drift of 
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the conduction carriers, the spontaneous sources remain essentially the 
same as in equilibrium. This will occur, for instance, when the sources 
originate primarily in the lattice, rather than in the system of conduction 
carriers, i.e, when the third term in (14) dominates over the imaginary part 
of the second term. Assuming that v is much smaller than the frequency of 
interest u> and comparing the two terms in (14), one arrives at the require- 
ment (e"/e' L ) ^> UpU/u 3 . Under this condition the correlation function of the 
spontaneous random sources is determined solely by the lattice and is given 
by [2,3]: 

= {jP(r)jf(T>))J( U -u/) (17) 

where it was assumed that the lattice, in the presence of the conduction cur- 
rent, remains close to equilibrium, at some temperature T. This assumption 
is quite common in the transport theory of metals and semiconductors. Thus, 
we arrive at a simple picture when the spontaneous fluctuations originate in 
the lattice while the conduction carriers and their drift only affect the subse- 
quent dynamics of those fluctuations, leading to modification of the spectral 
density of various fluctuating physical quantities. 

In what follows we shall be interested in the rotationless part of the 
fluctuating field. It is this part that determines the evanescent field close to 
the sample surface, as well as the short range correlations in the bulk of the 
sample. The rotationless electric field, E(r, t) = — V$(r, £), satisfies Eq. (3) 
with added random sources: 

On the right hand side of (18) appears the spontaneous random charge den- 
sity, p^, which is related to by the continuity equation, so that from (17) 
one obtains the expression 

<,<*V)^V)>^e>th|^(r-r<) (19) 

for the spectral density. In an infinite medium (18) can be Fourier trans- 
formed to obtain 

k 2 e{u, k)$(w, k) = 4ttp {s \uj, k) (20) 
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where e(u,k) is given by (12) or (13), depending on whether one keeps or 
not the thermal pressure term. 

Eq. (20), supplemented by the Fourier transformed Eq. (19) 

(pW(k)pW*(k')> w = ^(2vr)V>th|^(k-k') (21) 

enables one a straightforward calculation of correlation functions of various 
fluctuating quantities in an infinite medium (see Sec. 4). For finite bodies 
analytical treatment, in the presence of drift, becomes difficult (see Sec. 5 
for a specific example). 



4 Fluctuations in an Infinite Medium 

In this Section we consider fluctuations in the bulk of the sample, far from 
the boundaries. In this case one can assume an infinite medium and use Eqs. 
(20), (21), from which it immediately follows that the spectral density for 
the electric potential fluctuations is 

<<l(k)<r (k% = ^^^M 3 coth|^(k - k'). (22) 

Multiplying this expression by the factor k a kp and returning to real space 
one obtains the spectral density of the fluctuating electric field: 

r rl 3 h h ho P M r - r ') 
(*.(r W >. = 2 K coth| / (23) 

Let us make a short digression to discuss the equilibrium case when v = 
and e(cu, k) = e eq (u). In this case there is no point in making the assump- 
tion that the spontaneous random sources originate mainly in the lattice. 
One can rely on the fluctuation-dissipation theorem and use the full equilib- 
rium dielectric constant. For the Drude model, Eq. (14), one arrives at the 
expression 

{E a {r)^^{r J) ^equilibrium - Mj^—j—^ %f J (2~7t)3 ^ ' 

where, again, the double prime denotes the imaginary part of e cq (w). The 
integral in (24) exhibits an ultraviolet divergence, indicating a singularity for 
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r' — > r. For instance, setting j3 — a and tracing over a gives 



(E(r)£*(r')) w , equilibrium = ^j^L coth^5(r - r'), (25) 

which coincides with the second term in Eq. (88.24) of [2] or in Eq. (20.26) of 
[3]. To remove the divergence in (24) or (25) one must go beyond the Drude 
model and introduce spatial dispersion. One source of spatial dispersion is 
thermal pressure which resists any steep change of carrier concentration, thus 
introducing an ultraviolet cutoff in all integrals over k. All this is discussed 
in detail in Ref. 3 (see also Exercise 3.12.7 in Ref. 4), where another source 
of spatial dispersion, due to carrier diffusion, is also mentioned. 

We now return to (23). Writing |e(w,k)| 2 = \e'(u,k)\ 2 + |e"(w,k)| 2 and 
assuming that the imaginary part, e", is small, one can see that the important 
contribution to the integral comes from values of k in the vicinity of k c which 
satisfies e'(u,k c ) = (the pole contribution). But the condition e'(u,k c ) = 
is just the dispersion relation for the longitudinal plasma waves, in the 
absence of dissipation. The relation between plasma waves and fluctuations 
of the rotationless field becomes particulary clear for the spectral densities 
in the reciprocal space, such as 

(E a E;U = 2he" L coth^ ^ R ^yp (26) 



which is the Fourier transform of (23). Let us see how it works out for e(oo } k) 
given in (12). We take v — )■ and use the conditions kRn « 1, k- v << ou. 
For frequencies u close to cu p Eq. (12) simplifies to 

e(co,k) = e' L [F(cu,k)+tr ] }, (27) 

with 

F(u,,k) = 2 U -\- k - Vo -k*R* D , r ] = e \. (28) 
Up e L 

One can identify in (23) the combination i]/(F 2 + i] 2 ) which, for small rj, can 
be replaced by 7r5(F) yielding 

(E a Ep) ulti = -nhujp c °th^ S(u - u p - k • v - ^uj p k 2 R 2 D ). (29) 
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Comparison with (16) makes it clear that contribution to the spectral density- 
comes only from those regions in k-space where plasma waves can propagate. 
Transforming (29) back into real space, one can calculate the corresponding 
correlation functions. 

To elucidate the effect of drift on fluctuations let us consider the corre- 
lation function of the x component of the field, i.e., a = (3 = x, and fix uj 
somewhat below Co p . Furthermore, we chose Vo in the direction of the x-axis 
and consider correlations along re-direction, taking y = z = y' = z' = 0. 
Fourier transform of (29) then yields 



(E x (x)E* x (x'))„ = Tr^-coth-y _ e *.<*-* >._| (30) 
x 5{u - u p - k x v - l~u p k 2 R 2 D ), 



where only the essential arguments (x, x') in E x have been retained. In 
equilibrium this expression is zero, consistent with the absence of propagating 
plasma waves for u < Q p . Writing k 2 = k 2 + q 2 , where q is the transverse 
wave vector, and performing integration over q results in 

(E,(*)M = iJpr / kl + J F(uli K) m», *.))■ (3D 

where F(u,k x ) is given by Eq. (28), with k y = k z = 0, and the step func- 
tion selects the appropriate interval of k x in the integration region. For the 
integrand of (31) to be different from zero the condition 



v > R D \j2up{oOp - w) (32) 

must be satisfied. This is a necessary condition for propagation of plasma 
waves below cu p . The step function in (31) selects an interval of negative k x 
(i.e., opposite to the direction v ) such that 

7 - vV - 25 < \k x \R D < 7 + vV - 25 (33) 

where 

7=^, ' = (34) 
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The expression (31) takes the form 



(E x (x)E* x (x% 



h 



y+y/-y 2 -2S 



8ne' L R% 



du exp 



lu- 



ll 



D 



"fU 



(35) 



7-^-25 



This example demonstrates that drift can strongly modify the fluctuation 
spectrum in the vicinity of u p and, in particular, can lead to emergence of 
fluctuations at frequencies where there were no fluctuations in equilibrium. 
In addition, due to the oscillating factor in the integrand of (35), there will be 
oscillations in the field correlation function. Let us note that it is desirable to 
keep Vq smaller than Co p Re>. Indeed, £j p Rd is of the order of the carrier thermal 
velocity v t- If exceeds Vt, then heating of the carriers becomes appreciable 
and the hydrodynamic equation (7) would require some modification. Thus, 
we shall assume 7 to be smaller than unity. In Fig. 1 we give an example of the 
spectral density in Eq. (35), as a function of x — x', for 7 = 0.3, 5 = 0.01. The 
real part of the ratio 87re^i%(£ x (#)£*(x0) w /footh§^ = /(^) is plotted 

as a function of 



Rt 




10 15 

(x-x')/R D 
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Figure 1: The real part of the normalized spectral density (E x (x)E*(x')) u as 
a function of the distance 



5 Fluctuations near the surface 

In this section we study fluctuations of the evanescent electric field which 
exists close to the surface of a sample, due to fluctuating charges and carriers 
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inside the sample. We consider the simplest geometry of a sample occupying 
half space (z < 0) , while the second half (z > 0) is vacuum. A dc current 
is flowing in the medium and the conduction carriers are drifting in the x- 
direction with velocity v (Fig. 2). 



Z 




Figure 2: The sample at temperature T, with drifting carriers, is separated 
from the vacuum (T = 0) by a sharp boundary (the z = plane). 

Let us emphasize that in this section, unlike the previous one, the entire 
system (sample + environment) is out of equilibrium already in the absence 
of drift (v = 0). It is assumed that the sample is in local equilibrium, 
at temperature T, whereas the environment is "cold" (T = 0). In this case, 
which is often assumed in the studies of the electromagnetic field fluctuations 
[3,5], the sample is the only source of fluctuations so that no radiation is im- 
pinging on the sample from outside. Moreover, the zero-point fluctuations, 
which exist also in the vacuum, cancel out in the process of the electric field 
measurement. The latter belongs to the class of "absorption measurements", 
because some amount of energy must be diverted into the measuring device 
(the probe). This implies that instead of the symmetrized correlation func- 
tion for the random sources (Eq. (17) or (19)), with its characteristic factor 
coth|^, one should use the normally ordered correlation function (see [5] and 
references therein). The latter is obtained from its symmetrized counterpart 
by replacing ^coth|^ with the Planck function [exp{^) — 1] _1 = II(c<;,T). 
This replacement amounts to disregarding the zero-point fluctuations and it 
will be used throughout this section [19]. 

Because of the absence of translational symmetry in the z-direction it is 
not possible anymore to Fourier transform Eq. (18) in that direction. This 
complicates the analytic treatment considerably. One simplification, though, 
is that unlike the case of the infinite medium no ultraviolet cutoff will be 
needed in the present geometry. Therefore we will discard thermal pressure 
altogether, and use the permittivity tensor (11) without the last term. Fur- 
thermore, we will keep the k-dependence due to drift only in the component 
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e xx , thus arriving at a diagonal permittivity tensor: 

,2 



1 - 



+ i- 



[u - k x v + iv) (u - k x v ) e 



(36) 



-yy 



1 - 



+ i- 



Eq. (18), Fourier transformed in the x, y-plane, assumes the form: 



eo(w) 



d € xx (uj, k x ) ,2 , 7 2 
7T^7 H / s K x + K y 



dz 2 



$(h x , k y ] z) = A-Kp {s) (k X) k y ] z), (37) 



where tilde indicates the in-plane Fourier transform. The solution of this 
equation can be written in terms of a Green's function satisfying the equation 



(38) 



Since the right-hand-side of (37) differs from zero only inside the medium, 
while the solution of interest is outside the medium, we need g(z, z'\ q) for 
z 1 < 0, z > 0. The corresponding expression can be found, for instance, in 
Ref. 20: 

9(z, z'] q) 



-q{z-z>) 



(eo + l)5 

and the solution of (37) is written as 



$(k x , k y ] z) = 4tt / dz'g(z, z'; q)p (s) {k x , k y ] z), 



with 



e xx {uj,k x ) 2 2 
e o(w) \ 



1/2 



(39) 



(40) 



(41) 



The correlation function for the random sources p^ follows directly from 



Eq. (19), with the aforementioned replacement of ^coth 
U(u,T): 

(p^(k x ,ky]Z 1 )p^(k' x ,k' y ]Z 2 )} ul = 



hi) 
2T 



by 



(42) 



he'[U(uj,T)5(k x -k' x )5(k y -k' y ) 



d 2 



(k 2 x + k 2 y )5( Zl - z 2 ) + q^q^Kzi - z *) 
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With the help of (40), (42) and transforming back to real space, one obtains 



f° 

/ dzi [(kl + kl)g(z,z 1 ;q)g*(z',z 1 ;q) 4 

J —CO 



dg(z,zi;q) dg*(z',zi;q) 
dzi dzi 



• (43) 



From this expression one can calculate correlation functions for various com- 
ponents of the electric field. We limit ourselves to the ^-component of the 
field and consider correlations in the x-direction, for fixed y, z. The resulting 
function depends only on x — x' and z, and is given by 



f" 
1 1 




(E x (x,y,z)E* x (x',y,z)) 0J = 4h L U(ou,T) x 

foI^J + J- 1 

dk x dky 2 ik x {x-x>) ( k l + k l 1 \ 1 -2g' 2 (u) 

(2tt) 2 Ke V l# J 2g' ' 1 ' 

where q 1 is the real part of the quantity defined in Eq. (41). 

In the absence of drift, when e xx = e (u>), we have q = ^Jk 2 £ + k 2 y and the 
integral in (44) can be computed with the result 



m^Kw, „,,»<?> = r,w^ a[ ^j^§^ (45) 

where X = \x — x'\. Actually, more generally, one should replace in this 
equation e" L by €q(u). This is because in the absence of drift there is no 
need to make the additional assumption that the main source of spontaneous 
fluctuations is the lattice, rather than the conduction carriers. Eq. (45) is in 
agreement with the corresponding expression in [21] (our definition of spectral 
density differs by a factor 2tt from that in [21]). Setting in Eq. (45) X = 0, 
one observes that the spectral power of field fluctuations, (£ ,2 (^)) w , increases 
as z~ 3 when the surface of the sample is approached [3] (the divergence for 
z — > will be eventually cut off by some mechanism of spatial dispersion). 

It should be also noted that the factor \e (ui) + 1| in (45) is due to the 
surface plasmon wave which appears at the frequency satisfying the relation 
Reeo(u) — — 1. At this frequency, and provided that dissipation is small, 
a strong peak appears in the power spectrum [5]. On the other hand, the 
bulk plasmon frequency, which corresponds to Ree^iuj) = (i.e., oj ~ Co p ), 
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does not play any special role in Eq. (45). The frequency u p does become 
important, however, in the presence of drift, as we show next. 

Let us first take a closer look at the ratio e xx /eo which appears in expres- 
sion (41) for q. For v = this ratio is unity. For v ^ 0, however, it can 
become large if the frequency u is close to the bulk plasma frequency cu p . 
Indeed, e (w) will be then close to zero while € xx (lo, k x ) can be much larger 
due to finite vq. Taking v — > and assuming 



< |1- -^| < 1, (46) 

e L ' lo 



we have, from (36): 



e xx (u,k x ) 2(k x /k ) - (k x /h 



W) ' 2 V [1 - (k x /k W [ ' 



\2 



where 



r, = l-% ** = -■ (48) 

LO V 

In order to estimate the integral in Eq. (44) one has to find the relevant values 
of k x , making the main contribution to the integral. The factor exp(— 2q'z) 
in the integrand limits the value of q' to q' < 1/z, which in turn results in an 
efficient cutoff for k x , if z is not too small. In order to make a quantitative 
estimate we assume that the relevant region of k x corresponds to \k x \ <C ko 
and check the consistency of this assumption later. For \k x \ <C ko Eq. (47) 
simplifies to 

# e ^ = l_I^. (49) 
e r) k 

Since \r)\ <C 1, there exists a broad region of k x , 

\v\ko < \K\ < h, (50) 

such that the second term in (49) dominates and q (see (41)) can be approx- 
imated by 



k 2 - — 
y 7]k 



1/2 

(51) 
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This expression, together with the aforementioned condition q' < 1/z, implies 
that the effective cutoff for \k y \ is of order 1/z, whereas the cutoff for \k x \ is 
of order (Irjlko/z 2 ) 1 ^ 3 . Substituting this value into (50), and returning to the 
physical quantities ou, u p and v , gives the condition on z which is necessary 
for the picture to be consistent: 



\ W ~ U P \l/2 



^y /2 « 2 «^^. (52) 

V \UJpJ \u-uj p \ 

Using the cutoffs for k x and k y , one can now estimate the integral in (44) 
with the following result, for x' = x: 



(£2( ^»^ 2fi M^ nKr) ^ (53) 

Comparison between (53) and the corresponding result in equilibrium, i.e., 
Eq. (45) for X = 0, reveals that under the conditions specified above drift 
has a profound effect on the power spectrum of field fluctuations. The 
ratio between (53) and the corresponding equilibrium quantity is of order 
| a; — u) p \z/v , which due to (52) is much smaller than unity. Thus, our calcu- 
lation predicts a dip at the power spectrum for frequencies close to the bulk 
plasmon frequency uj p . 

The above estimate has been made under the condition given in Eq. (46), 
i.e., oj in Eq. (53) cannot be too close to u p . In order to approach the 
immediate vicinity of the bulk plasmon frequency one should replace the 
first inequality in (46) by the opposite one, \ui — Cj p \ <C u p {e" L /e' L ). Let us 
consider the extreme case u = u p . For this case Eq. (36), with v — > and 
\k x \v < u p gives 

6„/eo = l + 24^. (54) 

Note that this time it is essential to keep the imaginary part of this ratio. 
Moreover, for the effect of drift to be significant, the absolute value of the 
imaginary part must be much larger than unity, i.e., 

^ « \k x \ « ^. (55) 

The quantity q, Eq. (41), is now given by 

1/2 



q 



K + 2ik 



y 1 x ii 



e T to. 



v 



(56) 
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and the effective cutoff for k x in the integral (44) is of order (e'^ojp/e^voz 2 ) . 
Consistency with (55) requires z (vQ/Cj p )(e" L / e' L ) l l 2 . The integral in (44) 
can now be estimated and, for x' = x, we obtain: 

(El(x, y, z)) u ~ 2h { f 2 1 n(a;, T) (57) 

As expected, this value matches expression (53) at frequency oj such that 
\uj — Up\ ~ Cj p e" L /e' L . Eq. (57) gives the minimum, at a; = u p , of the dip in 
the spectral power. 

Our consideration has been limited to the case when the main contribu- 
tion to the spectral power, i.e., to the integral in (44), comes from k x which 
satisfy the condition \k x \v <C u>. This condition imposes a restriction on z, 
namely the first inequality in (52) (or its counterpart for frequencies very 
close to Cj p ). For z very close to the surface the integral in (44) is dominated 
by large k x , so that the condition |fc x |?;o <C w is violated. Considering values 
of k x larger than u p /vq while still neglecting the thermal pressure requires 
the condition v > Vt which is better to be avoided. 



6 Conclusion 

We have demonstrated the effect of carrier drift on thermal fluctuations of the 
electric field. Only the rotationless part of the field has been considered. This 
part dominates the short range correlations in the bulk of the sample, as well 
as the fluctuations in the near field sufficiently close to the sample surface. 
It has been shown that drift can significantly affect the magnitude and the 
correlation properties of the field fluctuations, especially at frequencies close 
to the bulk plasma frequency. Our main goal was to discuss the effect of 
drift in the simplest possible situation, making a number of idealizations 
such as the limit of collisionless semiconductor plasma {y — > 0) or keeping 
(in Section 5) the k-dependence due to drift only in the component e xx of the 
permittivity tensor. The latter simplification allowed us to avoid the intricate 
problem of the additional boundary conditions which might be needed in the 
more general case (see the book of Agranovich and Ginzburg in [17]). 

This work is in some sense complementary to [16] where a collision domi- 
nated transport regime was considered, i.e., it was assumed that the scatter- 
ing rate v is much larger than the frequency oo. In this regime the linearized 
Ohmic current is simply j = ctqE + env , where <r = e 2 n /mz/ is the dc 
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Drude conductivity. Adding to this the diffusion current —DVn, D being 
the diffusion coefficient, and eliminating n with the help of the continuity 
equation, one obtains instead of (13): 



e(u, k) =e' L 1 + 




(58) 



u — k • v + iDk 2 



where tm = e' L /4ir<7o is the Maxwell relaxation time. Neglecting e'[, one 
recovers from e(w,k) = the well known space charge waves, with the 
dispersion relation u = k • Vo — ^ — iDk 2 . These waves can strongly 
influence current fluctuations in semiconductors [12] and their impedance 
[22]. The surface counterpart of such travelling waves can be excited at 
the semiconductor-vacuum boundary [14] and can have a significant effect 
on the electromagnetic field fluctuations near the surface [16]. It would be 
worthwhile to further study the effect of these waves, under more general 
conditions that those assumed in [16]. 

One of the assumptions in the present work was that the spontaneous ran- 
dom sources of the fluctuations were not affected by the drift of the carriers. 
This is trivially so in the limit of collisionless plasma when the lattice remains 
the only source of spontaneous fluctuations. The assumption, though, can 
break down under more realistic conditions. It should be possible to relax 
this assumption. Indeed, for the current density fluctuations in the bulk, 
there is a well established theory [12, 13] for the case when the carriers are 
way out of equilibrium ("hot electrons") and the corresponding spontaneous 
random sources undergo a profound change. To our knowledge, so far there 
is no extension of the theory to the case of surface waves and their influence 
on the field fluctuations outside the sample. 

It is clear that drift, via its effect on the electromagnetic field fluctuations 
close to the sample surface, will influence also the Casimir-Lifshitz forces [23], 
as well as other related phenomena (heat transfer, noncontact friction). Cal- 
culation of the forces involves integration over frequencies, so that significant 
effect of drift can be expected only if the main contribution to the integral 
comes from the interval of frequencies in which the influence of drift on field 
fluctuations is strong. The calculation of the Casimir-Lifshitz forces in the 
presence of drift is beyond the scope of this paper. 

Finally, let us emphasize that the effects considered in this work are ex- 
pected to occur only in materials with relatively low carrier density (semi- 
conductors, ionic conductors or other types of "bad conductors"), where 
significant drift velocities can be achieved without destroying the sample. 
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